arXiv: 1509.02969v2 [astro-ph.HE] 23 Feb 2016 


Draft version February 24, 2016 

Preprint typeset using style emulateapj v. 5/2/11 


A TEST OF THE NATURE OF THE EE K LINE IN THE NEUTRON STAR LOW-MASS X-RAY BINARY 

SERPENS X-1 

Chia-Ying Chiang\ Edward M. Cackett^, Jon M. Miller^, Didier BarreU’'‘, Andy C. Fabian®, Antonino D’AP, 
Mighael L. Parker®, Sudip Bhattacharyya'’, Luciano BurderP, Tiziana Di Salvo^, Elise Egron^®, 

Jeroen Homan“, Rosario Iaria®, Dacheng Lin^®, and M. Coleman Miller^® 

^Department of Physics and Astronomy, Wayne State University, 666 W. Hancock, Detroit, MI 48202, USA 
^Department of Astronomy, The University of Michigan, 500 Church Street, Ann Arbor, MI48109-1046, USA 
®Universite de Toulouse, UPS-OMP, Toulouse, France 

^CNRS, Institut de Recherche en Astrophysique et Planetologie, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse cedex 4, France 
®Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, UK 
®INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica di Palermo, via U. La Malfa 153, 90146 Palermo, Italy 
^Department of Astronomy and Astrophysics, Tata Institute of Fundamental Research, Mumbai 400005, India 
®Dipartimento di Fisica, Universita degli Studi di Cagliari, SP Monserrato-Sestu, KM 0.7, 1-09042 Monserrato, Italy 
®Dipartimento di Fisica e Chimica, Universita di Palermo, via Archirafi 36, 1-90123 Palermo, Italy 
^®INAF-Osservatorio Astronomico di Cagliari, via della Scienza 5, 09047 Selargius (CA), Italy 
®®MIT Kavli Institute for Astrophysics and Space Research, 77 Massachusetts Avenue 37-582D, Cambridge, MA 02139, USA 
®^Space Science Center, University of New Hampshire, Durham, NH 03824, USA and 
^^Department of Astronomy and Joint Space-Science Institute, University of Maryland, College Park, MD 20742-2421 USA 

Draft version February 24, 2016 


ABSTRACT 

Broad Fe K emission lines have been widely observed in the X-ray spectra of black hole systems, and 
in neutron star systems as well. The intrinsically narrow Fe K fluorescent line is generally believed 
to be part of the reflection spectrum originating in an illuminated accretion disk, and broadened by 
strong relativistic effects. However, the nature of the lines in neutron star LMXBs has been under 
debate. We therefore obtained the longest, high-resolution X-ray spectrum of a neutron star LMXB 
to date with a 300 ks Chandra HETGS observation of Serpens X-1. The observation was taken under 
the “continuous clocking” mode and thus free of photon pile-up effects. We carry out a systematic 
analysis and find that the blurred reflection model fits the Fe line of Serpens X-1 significantly better 
than a broad Gaussian component does, implying that the relativistic reflection scenario is much 
preferred. Chandra HETGS also provides highest spectral resolution view of the Fe K region and we 
find no strong evidence for additional narrow lines. 

Subject headings: neutron star 


1. INTRODUCTION 


19951; iFabian et al.l 120091: iBrenneman et al. 

1201111. and 
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nass X-ray binaries (LMXBs) including stellar- 
black holes (iDone & Zvckilll999l: IMiller et al.ll2002l: 

Milled 120071: IReis et al. 1200911 and neutron stars (NS; 

Asai et al.l 120001: IBarret et al.l 120001: lOosterbroek et al. 
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:IDi Salvo et al.ll2005: Bhattacharvva & Strohmaver 
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An accretion disk is believed to orbit the cen¬ 
tral object, and a hard X-ray source, either a powerlaw 
continuum or a blackbody component (potentially the 
“boundary layer”), emits hard X-rays that illuminate 
the accretion disk. Atomic transitions take place after 
the high-energy photons are absorbed, resulting in a re¬ 
flection spectrum including several narrow emission lines 
and a broad feature peaked a round 20-30 keV which is 
known as “Gompton hiimn” (iLTsrhtman fc White! 119881 
George fc Fabianl ll991l: iR.oss fc Fabianlll993 : Matt et al 


s 


19931 IB OSS fc Eabia nI ]2nni~lG arcia fc KallmanI 12010 

Ballantvne et al.l 1201211 . 

The Fe Ka fluorescent line is the most prominent fea¬ 
ture in the reflection spectrum. When appearing in the 
X-ray spectra of AGN and LMXBs, the intrinsically nar¬ 
row Fe lines sometimes show broad, asymmetric pro¬ 


files which are generally believed to be shaped by a 
series of relati vistic effects induced from strong gravi¬ 
tational fields (jFabian et al.l[T9^ l2000l) . As relativis¬ 
tic effects are stronger in the area closer to the com¬ 
pact object, the line profile is sensitive to the inner ra¬ 
dius of the accretion disk. If the accretion disk extends 
to the innermost stable circular orbit (ISGO), one can, 
under certain assumptions, obtain an estimate of the 
black hole spin by measuring the inner radius of the 
accretion disk (jBardeen et al.l 1197^ . The Fe Ka line 
profile has been a powerful tool to measure black hole 
spin (e.g., iMiller et al.l 120021: iRevnolds fc NowakI 120031: 

[Brenn eman &: Reynolds! I2006t iMilleil l2007t iReis et al.l 

1200817120131 . The inner radius of the accretion disk in 
a NS system can be determined by the same method. 
The accretion disk in a NS system could be truncated 
by the stellar surface or the boundary layer between the 
disk and the NS, if the NS is larger than its ISGO, or 
by a strong stellar magnetic field. An upper limit of the 
stellar radius of a NS can be given by measuring the 
inner radius of the di sk, and hence help understand its 
equat i on of state (e.g. iPiraino et al.ll2000l : iGackett et al.l 
1200811 . iBhattacharvval (120111 1 reported more detailed cal¬ 
culations to show how future instruments can directly 
constrain NS equation of state models using relativistic 
disk lines. 

In NS systems, the Fe Ka lines are usually not as 
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prominent (EW ~ 100 eV) as those seen in AGN and 
black hole binaries (BHBs) due to extra continuum emis¬ 
sion from the boundary layer. It has been widely ac- 
cepted that relativis t ic Fe K lines are co mmon i n AGN 
(|Revnolds fc NowaS 1200311 and BHBs (IMilled 1200711 . 
though so me suggest line profiles to b e caused by warm 


absorbers (llnoue & Matsumoto 

12003(1 or Comptoniza- 

tion ((Laurent & Titarchukll2007 

). Nonetheless, whether 


Fe K lines in NS sys tems are r elativ istically broadened is 
still under debate. iNg et all (j20inll analyzed a number 
of XMM-Newton NS spectra and concluded that statis¬ 
tical evidence of asymmetric iron line profiles is lacking 
and the lines are broadened by Compton scattering in a 
disk co rona (iMisra fc Keinbhavi 119981 iMisra fc Sutarial 
1199911 . [Reynolds fc WilmsI (|200C 1 showed that the con¬ 
tinuum source required in the Compton scattering model 
to produce the broad iron line violates the blackbody 
limit. Although the calculation was based on the AGN 
case, small radii of the Compton clouds are still re¬ 
quired to maintain the high ionization level in NS sys- 
tems, which would make gravitational effects dominant 
(|Fabian et al.lll99^. Fu r therm ore, the pile-up correc¬ 
tion applied in iNg et'HI (j201flf l reduced the signal-to- 
noise ratio of the data and made it difficult to detec t 
relativistic iro n lines (see Figure 2 in iMiller et aLll^lOH . 
iCackett et al.l (1201011 examined a large sample of Suzaku 
NS spectra and found that the relativistic reflection sce¬ 
nario is much preferred. The iron line detections with 
Suzaku are less affected by photon pile-up effects than 
those of XMM-Newton, thus the conclusion that iron 
lines are asymmetric is likely more robust. A study of the 
effects of pile-up in X-ray CCD detectors bv iMiller et al.l 
(|2010ll showed that while pile-up can distort the Fe K line 
profi les, it tends to artificially narrow them (in contrast 
with INg et al.l 12010( 1 . 


Detections made with spectrometers that suffer no 
photon pile-up effects become important to determine 
the iron line profiles. The re-analyses of archival Bep- 
poSAX data of 4U 1 705-44 implied the existence of asym¬ 
metric Fe K l ines (iPiraino et a pj2003 iLin et all 120101 : 
ICackett et al.l 120121 : lEgron et ah 12013(1. NuSTAR sees 
an asymmetric line in Serpens X-1 ( Miller et al.ll20l'^ 
with line p roperties cons i stent with Suzak u measure¬ 
ments from ICackett et al.l (|2008L I2010L1201211 . The pile- 
up free detectors BeppoSAX and NuSTAR both reach 
the conclusion of relativistic iron lines. These instru¬ 
ments are, however, not capable of detecting possible 
narrow line components on top of the broad Fe K lines. 
Chandra HETGS is so far the only instrument that of¬ 
fers a pile-up free observation at high spectral resolution 
(approximate AE~30 eV at 6 keV). 

The neutron st ar LXMB Serpens X -1 was dis¬ 
covered in 1965 (jFriedman et al.l 11967( 1. Being a 
persistent, bright X-ray source. Serpens X-1 has been 
observed with majo r X-ray missions, includin g Einstein 
(IVrtilek et al.l[l9^ ASCA (lAsai et al.ll2nnnll. EXOSAT 
( Seon fc Mini 1200^ . BeppoSAX (lOosterbroek et al.l 
12001(1 . n^EGRAL (iMasetti et alJ I2004IV XMM- 
Newton iBhatte c harwa fc Strohmaved 12007( 1. Suzaku 
( Cackett et alll^OOM 12010(1 and recently with NuSTAR 
I Miller et al.l 1201^ It was als o detected in optica l 
( Hynes et al.l 12004(1 and radio (iMigliari et al.l 12004(1 
band. Relativistic Fe K lines have been reported 
several times in previous Suzaku, XMM-Newton and 


NuSTAR observations (i Bhattac harvva fc Stroh maved 
[20071: ICaciffitt et al.l [20Ml2010l: IMilleret all 1201, C h 
this paper we study the latest high-resolution Chandra 
HETGS observation, which is the longest Chandra 
grating observation of a neutron star LMXB to date. 
We present detailed data analysis and results in the 
following sections. The Galactic abs orption column Ah 
is as sumed to be 4.4 x 10^^ cm~^(lDickev fc LockmapI 
11990(1 with “wilm” abundances ([Wilms et al.l 1200^ 
throughout all our anal yses, which w ere done using the 
XSPEC 12.8.2 package ( Arnaudll 1996(1 . All errors quoted 
in the paper are given at the 90 per cent confidence 
level. 


2. DATA REDUCTION 

Serpens X-1 was observed with the Chandra High En¬ 
ergy Transmission Grating Spectrometer (HETGS) dur¬ 
ing 2014 June 27 — 29 and 2014 August 25 — 26 (Obs. 
ID: 16208, 16209), totaling a good exposure of ^ 300 
ks. The observation was taken using the “continuous 
clocking (CC)” mode. We reduced the data following the 
standard procedures using the latest CIAO V4.6 software 
package. The GG mode provides 2.85 ms time resolution, 
the observation is hence clear of photon pile-up effects 
and has negligible backgrounds. In this work we con¬ 
centrate on the HEG (high energy grating) data which 
covers the Fe line energy band. We found the HEG -(-1 
and —1 spectra to differ at the 5 — 10% level, especially 
around the area of a chip gap in the -1-1 spectrum. Thus, 
we tested a number of methods to improve the -|-1 spec¬ 
trum (see Appendix). However, the discrepancy cannot 
be completely eliminated, and thus we only use the HEG 
— 1 spectrum in this work. 

The spectra of each observation from June and Au¬ 
gust 2014 are similar and we combined them to form a 
long spectrum. All spectral bins from 2-8 keV have more 
than 30 counts per bin and thus no rebinning is required. 
We found a number of wiggles in the spectrum below 2 
keV that could not be modeled. Their locations matched 
those where there are significant sharp changes in the ef¬ 
fective area. We therefore use the 2.0 — 8.0 keV energy 
band in the following analysis. A restricted energy band 
is often used for the data taken under th e GG mode to 
avoid artificial instrumental artifacts (e.g. .lGacket t et al.l 
I2009bl: (Miller et alJI2011L (Ml IPe^aar et alJI20ilb 

The upper panel of Fig. [T] displays the 0.4 — 10.0 keV 
light curves of Serpens X-1 during the observation. It can 
be seen that there is a weak type I X-ray burst, which 
originates from t he thermonuclear burning of matter on 
the NS surfaces (I Wooslev fc Taaml ll976l: ILamb &: LambI 
119781: IStrohmaver fc BildstenI 120061^ in the later obser¬ 
vation (shown in blue data points). The X-ray burst 
lasted for a few hundreds of seconds and contributed 
^ 4.7 X 10"^ counts, which is only ~ 0.5 per cent of the 
total ~ 8.9 X 10® counts of the entire observation. Since 
the burst is such a small fraction of the total counts and 
casts no effects on the spectrum, we did not exclude the 
data during the X-ray burst. 

3. DATA ANALYSIS 
3.1. Continuum 

Gontinuum models for NS LMXBs can be degenerate, 
at least over narrow wavelength ranges, as we have here. 
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Table 1 

Best-fitting parameters for different continuum models. Model la/lb includes both DISKBB and BBODY 
components; Model 2a/2b composes of DISKBB and the powerlaw continuum; Model 3a/3b comprises of BBODY 
and powerlaw components. Both Model 3a and Model 3b gave most reasonable fitting parameters. 


Component 

Parameter 

Model la 

Model lb 

Model 2a 

Model 2b 

Model 3a 

Model 3b 

TBABS 

DISKBB 

BBODY 

Wh (10^^ cm“^) 
fc^disk (keV) 

^disk 

fcTbb (keV) 

Nbh (10-2) 

(0.44) 

80.ol®;® 

9 qQ”l"0-33 
^•'^^-0.19 
2.6 zb 0.1 

(0.44) 
1.42 ± 0.04 

es.sl®;! 

^•^^-0.44 

9 c+O-S 

(0.44) 
1.42 ± 0.01 
60.112-3 

(0.44) 

61-7t2.8 

(0.44) 

0.90lp;2° 
2.3 ±0.1 

(0.44) 

0 

i^-y^_o.oi 
2.4 zb 0.1 

POWERLAW 

r 

-Ypow 



0.94tO-23 
0.10 ±0.04 

0-68t®;2® 

U.UO Q Q3 

1.78 ±0.02 

0.88 ±0.03 

1.80 ± 0.02 
0.871^2 

GAUSSIAN 

Eline (keV) 

6.627®;®^ 

6.45_q Qg 

o.oz_0.07 

Q3 

fi 43+O Ub 

6.40+®-'-’2 


f7 (keV) 

Wga,u (10-3) 

EW (eV) 

-*-•^-0.5 

48+^^ 

0.61-0.11 

4 

^■°-1.4 

11^+42 

44 O 40 

0.41-0.03 

q 1 4-0.3 
*^•^-0.4 

yg-l-20 

0.61-0.03 
5.0 zb 0.6 

190+75 
-L4iO 29 

0 .41-0 .01 

q +0.3 
'^•-^-0.4 

74+16 

'^-14 

0.61-0.02 
5.3 ±0.6 
128111 

xVd.o.f. 


2052/1852 

2051/1852 

2072/1852 

2056/1852 

2184/1852 

2148/1852 



Figure 1. The upper panel shows the 0.4-10 keV light curve of 
Serpens X-1. The light curve was extracted from the first order 
HEG data. The gap between the observations has been discarded. 
Black points stand for data extracted from the first half of the 
observation, and blue ones for the second. There is a weak type 
I X-ray burst containing ~ 0.5 per cent of the total counts in the 
second half of the observation. The lower panel shows the hardness 
ratio, which is the hard (2-10 keV) light curve divided by the soft 
(0.4-2.0 keV) light curve. 

They can be equally well explained by different contin¬ 
uum models consisting of a disk blackbody, a blackbody¬ 
like component (to fit the boundary layer e mission) and a 
powerlaw or Co mptonized component (e.g., lBarretl[2001l : 
iLin et ani2007D . We test three different continuum mod¬ 
els, which are a disk blackbody plus a blackbody (Model 
1), a disk blackbody plus a powerlaw (Model 2), and a 
blackbody plus a powerlaw (Model 3), in order to find 
out necessary spectral components required to interpret 
the spectrum of Serpens X-1 in this observation. As 
mentioned in section [21 we use the HEG — 1 over the 
2.0 — 8.0 keV energy band. The Galactic absorption is 
modelled by TBABS in XSPEC, the accretion disk black¬ 
body emission by DISKBB, and the thermal emission 
from the NS boundary layer by BBODY. 

An iron emission line has been clearly detected in the 
HEG —1 spectrum, and we start by modeling the fea¬ 
ture with a Gaussian line (GAUSSIAN in XSPEG). In 
all our fits, we restrict the line energy to be in the range 


of 6.4 — 6.97 keV, where neutral/ionized iron Ka lines 
can only appear. However, given the relatively narrow 
energy range of the continuum, we find that the width 
(cr) of the Gaussian tends to be large values (1 - 1.5 
keV), and large normalizations giving equivalent widths 
(EWs) greater than 500 eV, even when fixing the energy 
of the line to 6.4, 6.7 or 6.97 keV. Since this is not seen 
in neutron star LMXB spectra, we restrict the width of 
the Gaussian based on previous fits to Serpens X-1 For 
instance, the broadest width that INg et all (I2010D get 
when fitting a Gaussian to AMMdata of Serpens X-1 is 
a = 0.27t°;lt keV. Fits to other archival spectra of Ser X- 
1 tend to give slightly broader Gaussians than this. For 
instance, averaging all the fits in IGackett et al.l (120121) we 
find an average Gaussian width of cr = 0.61 keV. Using 
these as a guide, we try spectral fits with (a) tr restricted 
to be less than the upper limit of 0.41 keV from Ng et 
ah' and (b) cr restri cted to be less than 0.61 keV from 
IGackett et al.l (I2012D . 

Including a broad Gaussian component in the iron line 
energy band significantly improves the fit (the smallest 
improvement was ~ 150 lower with three fewer 

d.o.f), confirming the clear detection of the broad iron 
emission line (also see differences between model residu¬ 
als in Fig. HD. In Tabled] we show the fitting results of 
the three possible continuum models. Note that we also 
tried fitting with all three continuum components, but 
found that the power-law index became unconstrained 
in those fits. 

Models lb & 2a-b give comparable quality fits, how¬ 
ever, all result in unphysical parameters. In Model 
lb, when a larger a is allowed than Model la, we get 
an unusually high blackbody temperature (see Table 
dD- Furthermore, when replacing the Gaussian with a 
DISKLINE, we get an unphysical blackbody tempera¬ 
ture (> 5 keV), which is insensitive to the data (the 
peak of the blackbody is well outside of the HEG en¬ 
ergy range). A typical blackbody temp erature seen in 
NS systems is generally below ^ 2.5 keVdLin et al.l 120101 

120121 : iPiraino et al.ll201^iDi Salvo et al.ll2015h . In Model 

2a-b the photon indices for the power-law are very hard 
(r = 0.94), which is also unexpected as photon indices of 
NS usually fall into the 1.5-3.0 range (IGackett et al.ll2010l . 

120121 : lEgron et al.ll2013l : IDi Salvo et al.ll2015D . Therefore, 

while Model 3a-b do not give the best fit statistically, we 
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use that model for the continuum in the following anal¬ 
ysis. 

It is interesting to note that both models la and lb 
fit equally well, yet have quite different Gaussian widths. 
This can be viewed as evidence that the broad line is 
asymmetric, since a relativistic line with a narrow core 
and broad red wing can be well fit by t wo Gaussians 
- one reaso nably narrow and one broad (iGackett et al.l 
I2008L l2012tl . Fits with two Gaussians to Suzaku spec- 
tra of Ser X-1 yield widths of 0.14 keV and 0.64 keV 
(jCackett et al.ll2?iT^ . comparable to the widths in mod¬ 
els la and lb. 

3.2. Iron Line 

An iron emission line was clearly detected in the spec¬ 
trum, implying that a reflection component is present 
in the system. The iron line of Serpens X-1 seems to 
extend for at least 1 keV, indicating possible signatures 
of relativistic effects similar to those seen in stellar-mass 
black holes and AGN. We test the nature of the line by 
fitting several different models. If the line is broadened 
by relativistic effects, it should be better fitted by a rela¬ 
tivistic line model rather than a broad Gaussian. Given 
the high spectral resolution of the Chandra HEG, we also 
have the opportunity to test whether there are any nar¬ 
row line components that contribute to the line shape 
that are otherwise unresolved with other detectors. 

First, we tested the relativistic reflection model 
by fitting the Fe e mission line using the DISKLINE 
(|Fabian et al.l 1198911 model in XSPEG. We replaced 
the broad Gaussian component in Model 3a/3b with 
DISKLINE to build a new model (hereafter Model 4a). 
See Table [2] for best-fitting values. This model fits signif¬ 
icantly better than a broad Gaussian, with an improve¬ 
ment of Ax^ = 133 for two more degrees of freedom. 
We hnd an inner disk radius of 7.7 ± 0.1 Rq (where 
i?G = GM/(?'), and an inclination of 24 ± 1 degrees. We 
show the iron line profile and the best-fitting DISKLINE 
model in Fig. |3l The DISKLINE component clearly 
fits the Fe line very well. Other relativistic line models 
such as LAOR and RELLINE were also tested, and all 
of them fit the Fe line as well as DISKLINE does, with 
similar parameters. Note that if we use Model la/lb for 
the continuum instead of Model 3a/3b, we get consis¬ 
tent DISKLINE parameters, and a large improvement in 
(Ay^ ~ 72 — 73), though because of the narrow en¬ 
ergy range we get a high blackbody temperature which 
is unconstrained. 

Next, we tested if the line could be fitted by two nar¬ 
row Gaussian components (Model 4b). The line width a 
of each narrow Gaussian component was set to be zero. 
Model 4b gives a significantly worse fit than Model 3b, 
indicating that narrow lines alone cannot fit the data. 
We next tested for the presence of narrow lines in ad¬ 
dition to a broad component (Model 4c), which gave a 
comparable fit than Model 3b, but a significantly worse 
fit than the relativistic line (Ay^ = 134 higher, with two 
fewer d.o.f.). A model including three narrow lines was 
tested but did not improve the fit. We also tried to set 
the line energies of the narrow Gaussian components to 
be those of the Fe XXV (6.67 keV) and Fe XXVI (6.97 
keV) lines, and again this did not yield a better fit. The 
equivalent widths of the narrow Gaussian lines in Model 
4b and 4c are all small (EW ~ 3 — 5 eV). We also tested 



2 3 5 8 

Energy (keV) 


Figure 2. The top panel shows the unfolded spectrum of Ser¬ 
pens X-1 with decomposed model components. We display 
the data/model ratios of the continuum model tbabs*(bbody + 
powerlaw), Model 3 and Model 3a in lower panels. 

using narrow Gaussian components with physical upper 
limits of the line widths as well. Assuming a narrow line 
originates from outer part of the accretion disk and is 
broadened by thermal effects in a ^ 10^ K gas, the line 
width should be ct < 0.007 keV. This again made no 
improvement in the fit. From the series of tests, we con¬ 
clude that there is no strong evidence of narrow emission 
lines. We present detailed 6.0-7.0 keV data/model ratio 
in the lower panel of Fig. |4l It can be clearly seen that 
Model 4a fits this band very well and no obvious emission 
or absorption features are present. 

In conclusion, of all the models we tried to fit the Fe 
K line in Serpens X-I, we find that a relativistic line 
model fits significantly better than any others, which in¬ 
dicates that the Fe line profile is caused by relativistic 
broadening and no narrow line components are required 
to explain the spectrum. The result acts to validate the 
relativistic reflection scenario. 

3.3. Relativistic Reflection Model 

As the Fe emission line is best interpreted by the rel¬ 
ativistic reflection scenario, we replace the DISKLINE 
component in Model 4a with a blurred broadband re¬ 
flection model. The reflection model self-consistently ac¬ 
counts for not only the Fe K line, but other lines and con¬ 
tinuum emission expected due reflection, and broadening 
due to Comptonization based on the ionization parame¬ 
ter. The reflection model is then blurred by relativistic 
effects. 

In NS systems, the accretion disk may be illuminated 
by the thermal emission coming from the boundary layer 
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Table 2 

Best-fitting spectral parameters, testing different models for the Fe K emission line. Model 4a 
tests a relativistic line, while models 4b tests narrow lines only and 4c tests a combination of a 
broad Gaussian and narrow lines. The relativistic line model is by far the best fit. Note that 
the definition of emissivity index /3 in the DISKLINE model is e(r) = (where r is the disk 
radius) instead of e(r) = and here we quote which is equal to —/3. 


Component 


Parameter 

Model 3b 

Model 4a 

Model 4b 

Model 4c 

TBABS 


A^h (10^^ cm“^) 

(0.44) 

(0.44) 

(0.44) 

(0.44) 

BBODY 


feTbb (keV) 

0.90l°;“ 

0.89l°;“i 

0.89 ±0.01 

0.90 ±0.01 



A'bb (10-2) 

2.4 ±0.1 

2.4 ±0.1 

2.2 ±0.1 

2.3 ±0.1 

POWERLAW 


r 

1.80 ±0.02 

1.74 ±0.02 

1.75 ±0.02 

1.78 ±0.01 



-Npow 

n 07+0.03 
0.02 

n +0-03 

U.oo o.o3 

0.88 ±0.01 

GAUSS 1 (broad) 


-^linel 

6.40+“'J^ 



g 4-I-U.U3 



(T (keV) 

0.61-0.02 



0.41-0 .01 



Afgau (10-3) 

5.3 ±0.6 



2.7 ±0.3 



EW (eV) 





GAUSS 2 (width= 

=0) 

-^line2 



6.57 ±0.01 

6.57 ±0.01 



Afgau (10-4) 



cy -| -1-0.6 

^•-*--0.4 

1.4 ±0.5 



EW (eV) 




stl 

GAUSS 3 (width= 

=0) 

-^lineS 



6.67 ±0.01 

6.67 ±0.01 



Afgau (10-4) 



1 O+0.5 

1.2 ±0.5 



EW (eV) 



c + l 
^-2 

3±1 

DISKLINE 


Eline (keV) 


6.97-0.02 





<? 


K gk+0.51 

O-OO_0 43 





Rin (RM/c 2 ) 


7.7±0.1 





inclination 


24 ± 1 





IVdiskline (lO-^) 


5.8 ±0.5 





EW (eV) 


149 ± 15 



/d.o.f. 



2148/1852 

2015/1850 

2286/1851 

2149/1848 



3 4 5 6 7 8 

Energy (keV) 


Figure 3. The figure shows the Fe line profile of Serpens X-1. The 
spectrum is modeled as tbabs(bbody + powerlaw + diskline), 
and the plot is produced by setting the normalization of DISKLINE 
to be 0. The plot is visually binned with “septplot rebin” in 
XSPEC. The bin sizes through the Fe line are ^ 0.042 keV. The 
red line shows model. It can been seen that the DISKLINE model 
fits the asymmetric line very well. 


of the NS or by a power-law continuum, resulting in re - 
fleeted emission (|Cackett et al.l[2M^ iD’Ai et alJ 2M^ . 
The continuum of Serpens X-1 in this observation is dom¬ 
inated by the po wer-law component, a nd we use the RE- 
FLIONX model (|Ross fc Fabianll2005li . which calculates 
the broadband reflection spectrum from the accretion 
disk illuminated by a power-law continuum. The con¬ 
volution model we use to account for relativistic effects 
is the KDBLUR kernel. We assumed the outer radius 
to be 400 Rq. The iron abundance Ape was set to vary 
in the range between 1 and 4 times of solar value. The 
model fits the spectrum well, and we report the best¬ 
fitting parameters in Table |3l The values of the inner 
radius i?in and inclination angle i we obtained are very 



Figure 4. The figure shows the 6.0-7.0 keV data/model ratio 
fitted using the DISKLINE model. The data are unbinned with 
bin sizes of 0.008 keV. It can be seen that no obvious emission 
or absorption features presenting in the figure. 

similar to those given by Model 4a. The model using RE- 
FLIONX yielded a better fit than Model 4a (Ay^ ^ 50 
lower with one fewer d.o.f.). 

We also tested the XILLVER reflectio n model 
(|Garcfa fc Kallmanl l2010bt iGarefa et al.l 1201311 . While 
the fit using XILLVER yielded a worse fit than that us¬ 
ing the REFLIONX grid, but the best-fitting parameters 
are comparable to that of Model 4a. In Table [3] it can be 
seen that models using different reflection grids gave con¬ 
sistent results. Parameters of KDBLUR obtained from 
both models are fairly similar, and both gave a inner 
radius of i?in ^ 7 — 8i?G and a low inclination angle 
{i - 30). 

4. DISCUSSION 

We analyzed a 300 ks G/iandro/HEG observation of the 
NS LMXB Serpens X-1. We fit a number of models to the 
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Table 3 

The table lists fitting results of the relativistic reflection models. It 
can be seen that REFLIONX and XILLVER yield similar parameters. 


Component 

Parameter 

REFLIONX 

XILLVER 

TBABS 

Afn 

(0.44) 

(0.44) 

BBODY 

fcTbb (keV) 

0.86 ± 0.01 

0.87 ±0.01 


iVbb (10-2) 

2.7 ±0.1 

2.6 ±0.1 

POWERLAW 

r 

1.68 ± 0.03 

1.75 ± 0.03 


-^pow 

0.56 ± 0.06 

n 'Tn + U.U2 
O.IO_o.04 

KDBLUR 

9 

. 1+1-8 
^•1-1.0 

4 4+1-1 


Rin (Rg) 

7 1+1.1 

^•^-0.2 


i (deg) 

29 ± 1 

33 ± 1 

REFLECTION 

Xpe 

l.00l°;5® 



« 

270t|° 



IVref 

G.Sfll X 10“® 

0 iq+0-64 
u.iy o.i5 

X^/d.o.f. 


1968/1849 

2012/1848 


2 — 8 keV HEG —1 spectrum to examine the nature of the 
Fe emission line, and find that the origin of the line is best 
explained by relativistically broadened reflection. Fitting 
broadband reflection models implies an inner radius of 

7 — 8 i?G and a low inclination of * 25° — 35°. In the 

following we discuss the robustness of the line parameters 
and compare our results with previous literature. 

4 . 1 . Choice of continuum model 

It is difficult to constrain the continuum using a re¬ 
stricted 2 — 8 keV energy band, and we choose the model 
with most reasonable fitting parameters (a blackbody 
and a powerlaw) to be the continuum in this work. In 
fact. Serpens X-1 has only been observed in the soft 
state, and the powerlaw component is usually weak. To 
explain the spectrum of a NS in the soft state, a con¬ 
tinuum composed of a disk blackbody component con¬ 
tributed by the accretion disk and a blackbody com¬ 
ponent possibly caused by the thermal emission from 
the boundary layer (Model la/Ib) is more likely. If re¬ 
placing the Gaussian component with the DISKLINE 
model in Model la/Ib and re-fitting the spectrum, we 
still obtain DISKLINE parameters similar to those of 
Model 4a (Elune = 6.93 ± 0.04 keV; emissivity index = 
-5.2t°j; i?in = 7.2 ± 0.1 Rq; i = 26 ± 1°). We also 
conduct the same test on Model 2a/2b, and find that the 
choice of continuum does not affect the parameters of the 
DISKLINE component. 

Assuming the continuum of Serpens X-1 is soft and 
dominated by the thermal emission from the boundary 
layer, the illuminating source is then the blackbody com¬ 
ponent. In order to further test the broadband rel¬ 
ativistic reflection model with a soft continuum (disk 
blackbody plus b lackbody), we use the BBREFL grid 
(|Ballantvnell2004 reflection calculated assuming a black¬ 
body component to illuminate the accretion disk) in¬ 
stead of REFLIONX to account for reflection. The model 
composed of a soft continuum and relativistic reflection 
(KDBLUR*BBREFL) gives the best-fitting inner radius 
i?in = 8 .1^5'2 Rq and inclination angle i = 34 ^3°, which 
are consistent with those obtained using a harder contin¬ 
uum with the REFLIONX/XILLVER grids. 

4 . 2 . No additional narrow line components 

In section 13.21 we test various models to fit the Fe K 
line and show that the line can simply be fitted by a rel- 



Figure 5. The figure shows the contour plot of the emissivity 
index against the inclination angle, which was calculated from the 
REFLIONX model. The contours are plotted at 67% (black), 90% 
(red) and 99% (green) confidence levels. It can be seen that there 
is no claer degeneracy between emissivity and inclination angle. 

ativistically broadened reflection component. Given the 
unique spectral resolution of Chandra HETGS, we have 
an opportunity to test for the presence of any narrow 
lines, in addition to the broad line. We find that includ¬ 
ing narrow lines in addition to a broad Gaussian gives 
a worse fit than the relativistic line alone. We include 
narrow Gaussian components with reasonable line ener¬ 
gies (6.4 keV, 6.67 keV and 6.97 keV) to Model 4a and 
the broadband relativistic reflection models to examine 
the existence of narrow lines. Narrow components at 6.4 
keV and 6.67 keV do not improve the fit of Model 4a, 
and the narrow Gaussian line at 6.97 keV improves the 
ht marginally (Ay^ ~ 3 lower than Model 4a). The nar¬ 
row lines at 6.67 keV and 6.97 keV improve the fit of the 
REFLIONX model marginally (Ay^ ~ 5 — 10 lower than 
the REFLIONX model), but equivalent widths of these 
lines are low (EW ^ 2 eV). Furthermore, examination of 
Figure |4] shows no evidence for narrow lines in the resid¬ 
uals of the relativistic line fit. Hence, we conclude that 
there is no strong evidence of narrow components. 

4 . 3 . Comparison with Previous Work 

iBhattacharvva fc Strohmaved (j2007D . iGackett et al.l 
(j20ir)l ) and iMiller et al.l (|2013|) reported the existence 
of relativistic Fe K lines in Serpens X-1 using XMM- 
Newton, Suzaku and NuSTAR dat a, respectively. In 
[Bhattacharvva &: Strohmaved (120071 ) the continuum was 
modeled as an absorbed Gomptonization (compTT) plus 
a disk blackbody, while in the later two pieces of work 
a full continuum (a disk blackbody, a blackbody plus 
a pow erlaw) was used. IBhattacharvva fc StrohmaveH 
(l2(^ used the LAOR model to fit the Fe K line, while 
IGackett et al.l (I2010D and IMiller et al.l (I2013D used the 
blurred BBREFL and R EFLIONX models to account 
for re lativistic reflection. Bhattacharvva fc StrohmaveH 
(|2007D obtained a inclination angle of f ~ 4 0°—50°, which 
is hig her than the best -htting valu e s of IGackett et all 
dMff) {i < 25°) and IMiller et al.l (I2013n (i < 20°) 
and this work {i ^ 30°)^ IGackett et al. (|20ICT also 
analyzed the same XMM-Newton data set used in 
[Bhattacharvva fc StrohmaveH (|2007f) and found that low 
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i nclination (z < 30°) is pref erred. _ 

iCackett et al'l (|2010f) and iMiller et al.l (I2013D both ob¬ 
tained a low emissivity index of g ^ 2.3, while a higher 
value of <7 4 — 5 is required in this work. We note 

that incli nation and emissivit y can be degenerate (see 
Fig. 8 in ICackett et al.l [201011 . but the degeneracy was 
not found in this work (see Fig. [SJ. The emissivity in¬ 
dex depends on the Fe K line profile, and it is likely that 
the line profile is slightly different from those of previ¬ 
ous observations because the better spectral resolution 
of Chandra grating data in this energy band. Regard¬ 
less, our fits suggest a small inner radius, consistent with 
previous findings. In this work, the broad iron line pro¬ 
file implies an inner radius of ^ 7 — 8 Rq. Previous 
analyses on XMM-Newton, Suzaku and NuSTAR data 
suggest the inner radius to be ^ 4 — 25 Rq, ~ 6 — 8 
i?G, and ^6 — 8.5 Rq, respectively. The XMM-Newton 
observation was taken under timing mode with a short 
exposure time, which causes some uncertainties in con¬ 
straining parameters. The later Suzaku, NuSTAR and 
current observations all give well-constrained, consistent 
measurements of the inner radius. Although different 
continuum models were used to fit the spectra of Ser¬ 
pens X-1 observed with various instruments at different 
times/fluxes, previous and current analyses all indicate 
the source to have low inclination angle and small inner 
radius. 

All previous observations of Serpens X-1 discussed 
above place the inner disk radius close to the ISCO. 
When comparing to the broader class of NS LMXBs, 
most other sources have an inner disk radius con¬ 
sistent with the ISCO, regardless of source lumi¬ 
nosit y (for instance see Figure 7 of ICackett et al.l 
uMa for a comparison over two orders of magni¬ 
tude in luminosity) . Clear exception s to this are 
SAX J1808.4—3658 (|Cackett et al.ll2009al: Paoitto et ahl 
HOOl), IGR J1 70480-2446 (IMiller et al.l I2011D . ^ 
GRO J1744—28 (|Degenaar et al.ll2014il . all of which are 
pulsars whose magnetic field must truncate the inner ac¬ 
cretion flow. In low/hard states, accretion flow models 
suggest that the disk should recess, be ing replaced by a 
hot radiatiyely inefficient flow (see e.g., iDone et 'mi2n07l 
for a reyiew). Howeyer, iron line studies of NS LXMBs 
in low states do not show large disk truncation. For 
instance, a NuSTAR obseryation of 4U 1608-52 at 1-2 
per cent of the Edd ington limit has a dis k that extends 
close to the ISCO (|Degenaar et al.l[2015D . while broad¬ 
band Suzaku obseryations of 4U 1705—44 at 3 per cent 
Eddin gton also shows the disk is not truncated at large 
radii l|Di Salvo et al.l 120151) . Thus, the majority of NS 
LMXBs, including Serpens X-1, show inner disk radii 
that are close to the ISCO, with no strong dependence 
on state or luminosity. 

5. CONCLUSION 

We analyze the latest long HETGS data of Serpens X-1 
and examine the nature of its Fe emission line. A ther¬ 
mal blackbody component possibly contributed by the 
boundary layer of the NS, and a power-law component 
provides a good fit to the continuum without unphysical 
parameters. By studying the Fe emission line in the spec¬ 
trum, we find that the relativistic reflection scenario is 
much preferre d , whic h is consistent with previous studies. 
ICackett et al.l (|2010D analyzed Suzaku data of Serpens 


X-1 and found relativistically-blu rred iron emission lines. 
The recent NuSTAR observation (|Miller et al.ll201^ con¬ 
firms the presence of the relativistic iron line, together 
with the Compton hump. In this work, we construct sev¬ 
eral models to test the relativistic reflection scenario and 
find that blurred reflection explains the Fe line profile 
significantly better than single/multiple Gaussian lines. 
Thanks to the remarkable resolving power of Chandra 
HETGS, the grating spectrum is capable of detecting 
narrow emission lines. In our analysis, no narrow line 
components are required, and must be weak if existent. 

The broad iron line profile implies a small inner radius, 
and we obtain an inner radius of ~ 7 — 8i?G- This sets an 
upper limit to the NS radius of ~ 15 — 17 km (assuming 
the mass of the NS is ^ 1.4 Mq). Eurthermore, a low 
inclination angle of ^ 25 — 35 degrees is found, which is 
consistent with the previous measurements. We also find 
that the choice of continuum does not affect the values of 
the line-related parameters, which further confirms the 
robustness of the fitting results. We conclude that the Ee 
emission line observed in the X-ray spectrum of Serpens 
X-1 is broad and shaped by relativistic effects. 

This work was greatly expedited thanks to the help 
of Jeremy Sanders in optimizing the various convolution 
models. CYC and EMC gratefully acknowledge support 
provided by NASA through Chandra Award Number 
GO4-15041X issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophys- 
ical Observatory for and on behalf of NASA under con¬ 
tract NAS8-03060. 
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APPENDIX 

DISCREPANCY BETWEEN THE +1 AND -1 SPECTRA 

We aim at analyzing data with the maximum signal-to-noise ratio, so seek to combine the HEG +1 and —1 spectra. 
The HEG ±1 spectra are, however, discrepant and not suitable for combination. In Eig. |6] we plot the detector 
effective areas and data/model ratio of the HEG ±1 spectra (+1 in black and —1 in red data points) fitted by a simple 
continuum tbabs* (diskbb + bbody + powerlaw). It can been seen that the ±I spectra disagree in most of the HEG 
energy band (at the 5 — 10% level), though both show similar Fe line profiles when the continuum is properly modeled. 
The wiggles and emission features shown in the —1 spectrum below 2.0 keV cannot be modeled and match changes in 
the effective area, thus are likely due to calibration uncertainties. There are obvious deviations over the 2 — 5 keV band 
between the ±1 spectra. It seems the HEG +1 spectrum suffers calibration problems, while the HEG —1 spectrum 
does not show unexpected features in the 2 — 5 keV energy band. The location of the zeroth order image on the chip 
determines at which energies the chip gaps lie. Here, a 1 arcmin y-offset was applied on the zeroth order in order to 
place it in a location to avoid any chip gaps near the Fe K band. This results in a chip gap at around 2.5 keV in 
the +1 spectrum. Use of such an offset is not widely reported previously, which may explain the issues here. The 
^^2.5 keV drop in the +1 spectrum matches with a large change in the effective area due to a chip gap (as marked by 
dashed lines in Fig. ii). Uncertain calibration around this gap clearly leads to the residuals. There also seems to be 
excess emission around the 3 — 4 keV. This has been a known issue for HETGS data in CC mode, probably causing 
by improper order sorting table (OSIP) or charge transfer inefficiency (CTI) corrections^. 

Different from the “time event (TE)” mode, in CC mode every CCD suffers from time-dependent CTI. The CTI 
correction relies on charge trap maps for each device to predict charge losses and correct them. Nevertheless, in CC 
mode charges are clocked continuously and possibly leads to time-dependent charge trap maps, and each CCD suffers 
from this effect. Inappropriate CTI correction may cause events to fall out of the OSIP. So far, alternate trap maps 
for CC mode are not found, but the Chandra calibration team provided a few possible methods to solve this issuqj- 
One can primarily use HEG —1 and MEG -1-1 orders only, or apply a custom OSIP to possibly fix the problem. 

We examined the order-sorting regions of each CCD by plotting the grating wavelength (wavelength times order) 
against the wavelength over the CCD wavelength (hc/ENERGY). Only HEG first order data have been used, and we 
show the order-sorting regions of chips S2 and S3, where the ^2 — 5 keV ±1 spectra were extracted from, in Fig. [7l 
Ideally data points should be evenly distributed along the y = 1 axis. Most events on chip S2 lie slightly above unity 

^ http://cxc.harvard.edu/cal/Acis/CaLprods/ccmode/ccmode_finaLdoc02.pdf 
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Figure 6. We plot the HEG +1 (in black) and —1 (in red) detector effective areas in the upper panel. The lower panel of the figure shows 
the Serpens X-1 HEG +1 (black data points) and —1 (red data points) data/model ratios. The continuum is modeled as tbabs* (diskbb + 
bbody + powerlaw). The plot reveals the poor calibration of the HEG +1 spectrum. The vertical dashed lines are intended to help guide 
the eye. 
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Figure 7. The order-sorting region of chips S2 and S3. The x-axis is the wavelength times order, and the y-axis the absolute value of 
x-value over the GCD wavelength. The figure shows the selected events during the order-sorting process. The mild curvature shown in S3 
might be the cause of discrepancy between the HEG ibl spectra. 


(and so as those on SI, S4 and S5), while the distribution of those on chip S3 shows mild curvatures, which might 
be the reason that the +1 spectrum is not consistent with the —1 spectrum. One possible way to improve the ±1 
spectral agreement is to modify the event file and bring the y-value in Fig. [3 close to unity. We first tried to correct 
the ENERGY column in the level=1.5 (after order-sorting) event file by dividing the column by the average y-value 
of each chip, and extract spectra from the corrected event file. We also tried a more sophisticated method, which is 
modifying the ENERGY column node by node, i.e., applying a spline fit to make every event lies on y = 1. However, 
neither of the methods made the spectra noticeably better. We then tried to modify the level=l (before order-sorting) 
event file and re-run tg_resolve_events to process order-sorting on the corrected event file, and extracted a new 
spectrum from the new event file after order-sorting. Modifying the level=l event file does give different spectra, but 
it seems it causes problems to order-sorting and the spectra are clearly not corrected properly. 

Another method to tackle this issue is to widen the order-sorting window. When data are taken under the CC mode, 
the Y position of an event is not known, and the order selection can be tricky. In some cases the spectra can be smoothed 
out by including events fell out of the OSIP. When running tgjresolve_events, by setting the script to disable the 
original OSIP file and indicating numbers of parameters “osort_hi” and “osort Jo”, one can customize the size of the 
order-sorting window. We set both osortjii and osortJo to be 0.2 (including events fell in the regime 0.8 < y < 1.2 
in Fig. [T]) to widen the event order-sorting window. We find that this does not help eliminate discrepancies between 
the ±1 spectra, but on the contrary, makes the issue worse. By widening the order-sorting window, more events are 
selected and hence spectra with higher fluxes are created. The 3-4 keV excess emission in the -1-1 spectrum turned 
to be larger than the original spectrum before correction. In this case, widening the order-sorting window does not 
effectively solve the problem. 

We look for a correction that would remove the 3-4 keV excess from the HEG -1-1 spectrum. As widening the 
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order-sorting would increase the flux of this energy band, narrowing the window might give a correction that we need. 
We And that setting osort_hi=0.2 and osortJo=0.04 (including events fell in the regime 0.96 < y < 1.2) improves 
the agreement of the ±1 spectra. The corrected -1-1 spectrum does not completely match the —1 spectrum, but the 
flux level between 2 — 5 keV is much more similar than the spectrum before correction. Although unfortunately it is 
still not good enough to combine the ±1 spectra and achieve maximum signal-to-noise ratio, we find that applying a 
SMEDGE component in XSPEC with negative optical depth to the -1-1 spectrum when fitting significantly improves 
data agreement over the 2 — 5 keV energy band. The SMEDGE component can mimic the sharp feature caused by 
the change in effective area and reduce the residuals. Yet it seems that different continuum models are required to fit 
the -1-1 and —1 spectra, thus we only present results of the —1 spectrum in the paper. Nevertheless, fitting the -1-1 
spectrum including the SMEDGE and allowing for a different continuum model than the —1 spectrum, results in the 
same conclusions and consistent parameters for the relativistic iron line. 



